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Abstract 



We examine hysteresis in cyclic voltammetry of Br electrosorption onto single-crystal Ag(100) 
by dynamic Monte Carlo simulations. At room temperature, this system displays a second-order 
phase transition from a low-coverage disordered phase to a doubly degenerate c(2 x 2) phase. 
The electrochemical potential is ramped back and forth across the phase transition, linearly in 
time, and the phase lag between the response of the adlayer and the potential is observed to 
depend on the sweep rate. The hysteresis in this system is caused by slow ordering/disordering 
kinetics and critical slowing-down. 

Hysteresis is a widely occurring phenomenon in which a system is driven too fast to remain near 
equilibrium, resulting in a phase lag between the response and the oscillating driving force. Some 
common examples include magnetic systems in oscillating external fields Jj], ffl, electrochemical ad- 
sorption in an oscillating electrochemical potential H |], ||, and nonlinear media under oscillating 
stress. As an example of hysteresis at a second-order phase transition, we here present prelimi- 
nary results of dynamic Monte Carlo simulations of the electrosorption of Br onto single-crystal 
Ag(100) under a time-varying electrochemical potential, known as cyclic voltammetry (CV) in the 
electrochemical literature. 

In equilibrium at room temperature, this system displays a second-order phase transition in 
the Ising universality class from a low-coverage disordered phase at more negative potentials to a 
doubly degenerate c(2 x 2) ordered phase at more positive potentials ||. As the system is driven 
back and forth across this phase transition, hysteresis is observed H). Unlike magnetic systems 
|0, |^] and some electrochemical systems |Q , in which hysteresis is usually associated with metastable 
decay near a first-order phase transition Jj], Q|, hysteresis in the Br/Ag(100) system is associated 
with kinetic limitations to the phase ordering/ disordering processes and critical slowing-down in the 
neighborhood of the second-order transition. 

The model and the dynamic Monte Carlo algorithm are discussed in detail in Ref. || and are only 
briefly summarized here. In-situ X-ray scattering indicates that the Br adlayer is commensurate with 
the square Ag(100) lattice, with Br adsorbing at the four-fold hollow sites between the Ag atoms J|]. 
This arrangement suggests a lattice-gas treatment of the adlayer. We used an L x L square lattice 
of adsorption sites with periodic boundary conditions. The grand-canonical lattice- gas Hamiltonian 



where i and j denote sites on the lattice, c; denotes the occupation of site i (0 for empty and 1 for 
occupied), 52j<j denotes a sum over all pairs of sites on the lattice, <pij denotes the lateral interaction 
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energy of the pair and p is the electrochemical potential. The sign convention is such that 

<j>ij < denotes a repulsive interaction, and p > favors adsorption. 

The interactions are quite adequately described with a nearest-neighbor excluded-volume inter- 
action (infinite repulsion) caused by the large ionic radius of Br, plus a long-range dipole-dipole 
repulsion jH[ |J. The lateral pair interactions are H 

!— oo r = 1 

2 3/2 <Wr 3 V2<r<5, (2) 
r > 5 

where r is the separation of an interacting Br pair within the adlayer, measured in units of the 
Ag(100) lattice spacing, and cj> nnn is the value of the repulsion between next-nearest neighbors. 
The repulsion was truncated to simplify computation. We found <p nnn = —26 meV/pair by fitting 
equilibrium simulation isotherms to experimental adsorption isotherms ||. 

To simulate the microscopic dynamics, we used a thermally activated stochastic barrier-hopping 
scheme, which include adsorption, desorption, and nearest- and next-nearest neighbor diffusion. 
These processes were represented by transitions between initial and final lattice-gas states, I and 
F, respectively, through intermediate transition states of higher energy. These intermediate states, 
T\, cannot be represented by lattice-gas configurations, and the subscript A indicates the process 
(adsorption/desorption, etc.) which relates / to F. The energies of these transition states were 
approximated as Ji|, ||] 

U Tx = (Ui + U F )/2 + A x , (3) 

where Ut x , Ui, and U f are the energies of the states T\, I, and F, respectively, and A> is the barrier 
associated with the process A. 

We approximated the probability R(F\I) of making a transition from I to F during one time 
step by an Arrhenius rate [||, ||] 

R(m . _ p . „« p (_ exp (4) 

where fes is Boltzmann's constant, T is the temperature, and v is a dimensionless attempt frequency 
which sets the overall time scale of the simulation. For all simulations discussed here, the temperature 
was 290 K, corresponding to k^T = 25 meV. The nearest-neighbor diffusion barrier was A nn = 100 
meV and the next-nearest neighbor diffusion barrier was A nnn = 200 meV, consistent with ab initio 
calculations |Q. The adsorption/desorption barrier was A a = 300 meV, and we used v = 1 [||. Time 
was measured in Monte Carlo steps per site (MCSS). 

The nearest-neighbor excluded-volume interaction gives rise to two degenerate c(2 x 2) ordered 
phases, each occupying one of two sublattices, A and B. The sublattice coverage, 8a, is the fraction 
of occupied sites on sublattice A, and analogously for Ob- There are two observables of interest: 
the total Br coverage, 8 = (8,4 + 8s)/2, and the staggered coverage 8s = &a — 8b, which is 
the c(2 x 2) order parameter. Figure 1 shows equilibrium 8 and |8s| isotherms obtained by Monte 
Carlo simulation using a single-site Metropolis algorithm The transition occurs at a critical 
electrochemical potential, fi c = 180 ± 5 meV/particle J|]. 

Although hysteresis studies in magnetic systems usually use sinusoidally varying fields, CV ex- 
periments mostly use linearly ramped electrochemical potentials. For the hysteresis (dynamic CV) 
simulations, the Br/Ag(100) system was first equilibrated in the low-coverage phase at fix <C n c - 
Then fi was ramped linearly in time up to fi2 3> fJ-c at a rate of p meV/MCSS. After reaching the 
potential was ramped back down to /xi at the same rate ||. Because of the long simulation times, 
only a single cycle was simulated for each value of p. 

Figure 2 shows hysteresis loops, 8 vs. p, which are analogous to magnetization vs. field plots 
commonly shown for hysteresis in magnets |Q, B. However, for CVs, dO/dt vs. p is usually shown, 
since dQ/dt is directly measurable as proportional to the electric current through the electrochemical 
cell, Fig. 3. This figure also shows a quasi-equilibrium CV, obtained by differentiating the equilibrium 
coverage isotherm in Fig. 1 with respect to p. In the limit of zero sweep rate, the dynamic CVs 
should have the same shape and p placement as the quasi-equilibrium CV for the same system size. 
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The asymmetry between the positive-going scans (positive dQ/dt) and the negative- going scans is 
caused by the difference in the ordering (positive scan) and disordering (negative scan) dynamics. 
In the limit of zero sweep rate, positive and negative scans would be symmetric. Figure 4 shows the 
height, (d0/di) max , and position, p, p , of the sharp peak in the positive-going scans vs. p. 

In conclusion, we have studied the hysteresis of the Br/Ag(100) adlayer by dynamic Monte Carlo 
simulation. We observed asymmetry between the positive and negative-going scans, which is ex- 
pected to vanish in the limit of zero sweep rate. We also observed a power-law relationship between 
p _1 (d6/dt) max and p with a measured effective exponent of w —0.08. Future work will include 
slower sweep rates to test these observed behaviors. 
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Figure 1: Equilibrium Monte Carlo isotherms for L=32. A second-order phase transition between 
a low-coverage disordered phase and a c(2 x 2) phase with 6=1/2 is observed at fi c - The isotherm 
was generated from 10,000 independent samples for each value of /U. 
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Figure 2: Hysteresis loops at various sweep rates p. All curves show vs. /i for L=256, except for 
the solid curve, p=0, which shows equilibrium for L=32, Fig. 1. 
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Figure 3: Simulated CVs for L=256 at various sweep rates p. All curves show dQ/dt vs. /i, except 
for the solid curve, p=0, which shows dO/dp at equilibrium for L=32. All CV curves, except p=0, 
have been normalized by the sweep rate. 
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Figure 4: CV peak position and height vs. sweep rate, (a) shows the positive-going peak height 
with errors much smaller than the symbol size, (b) shows the peak displacement from /Lt c , fi p — fi c , 
and the distance in \x between 9=1/4 in the positive-going and negative-going scans in Fig. 2. 
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